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We extend our investigation of the normal modes of small adiabatic oscilla 



6 ■ 

tions of relativistic barotropic thin accretion disks to the inertial-pressure (p) 
modes. We focus here on the lowest frequency fundamental p-modes, those with 
no axial or vertical nodes in their distribution. Through a variety of analyses, 



we obtain closed-form expressions for the eigenfrequencies and eigenfunctions. 
These depend upon the luminosity and viscosity parameter of the disk, as well 
as the mass and angular momentum of the black hole, via detailed formulae for 
the speed of sound. The effect of a torque on the inner edge of the disk is also 
included. We compare the p-mode properties to those of the g- and c-modes. 
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Introduction 



All perturbations of equilibrium models of accretion disks can be described in terms 
of their normal modes of oscillation, which our group has been studying d uring the last 
decade, following the pioneering work of iKato k, Fukud (119801 ) [see iKatol (j200ll ) for a review]. 
Approaches involving instead more local (e.g., dispersion relation) considerations have been 
extensively used to study the propagation of waves in disks , with th e inclusion of buoyancy, 
atmospheres, and magnetic fields [see, e.g., Lubow fc Ogilviel ( 1998 ) and references therein]. 
Their classification of modes is similar to those for stars, and different than ours. Our focus is 
rather on global modes, and in the determination of eigenfrequencies and their applications. 
All the modes we study have discrete, real eigenfrequencies. Viscous damping (or growth) 
is assumed to introduce a relatively small imaginary contribution. 

Effects of general relativity trap most of these modes near the inner edge of accretion 
disks around compact objects. The strong gravitational fields that are required can be 
produced by a black hole or by neutron stars that are sufficiently compact (with a soft 
equation of state) and weakly magnetized to produce a gap between the surface of the star 
and the innermost stable orbit of the accretion disk. Although we shall not explicitly consider 
such neutron stars here, the results obtained will also apply to them to first order in the 
dimensionless angular momentum parameter a = cJ/GM 2 , since their exterior metric is 
identical to that of a black hole to that order. 



The subj e ct of 'relativistic diskoseismology' has been reviewed by IKato et al.l (119981 ) 
and IWagonerl (119991 ). In this paper we shall focus on the fundamental low frequency 
(LF) pressure-driven p-modes: those axisymmetric modes with no nodes in their ver- 
tical distribution and with the smallest possibl e eigenfrequency. We shall briefly com- 



pare 



these modes with the 'corrugat ion' (c) modes (ISi 



20011 ) and the 'gravity' (g) modes (IPerez et al.l Il997l ). An investigation of the other types 



bergleit. Wagoner fc Ortega- Rodriguez 



of p-modes, including those which can coyer th e entire disk, will be published separately 
( lOrtega- Rodriguez. Silbergleit fc Wagonerll200ll ). The study of p-modes will complete the 
last phase of our analysis of the normal modes of oscillation of 'standard' relativistic thin 
accretion disks. 
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2. Basic Assumptions and Equations 



2.1. Structure of the unperturbed accretion disk 

We take c = 1, and express all distances in units of GM/c 2 and all frequencies in units of 
c 3 /GM (where M is the mass of the central body) unless otherwise indicated. We employ the 
Kerr metric to study a thin accretion disk, neglecting its self-gravity. The stationary (d/dt = 
0), symmetric about the midplane z = 0, and axially symmetric (d/ dtp=0) equilibrium disk is 



taken to be described by the standard relativistic thin disk model (INovikov fc Thorndll973 
Page fc Thorndll974j ). The velocity components v r = v z = 0, and the disk semi-thickness 
h(r) ~ c s /Q <C r, where c s (r, z) is the speed of sound. The key frequencies, associated with 
free-particle orbits, are 



ft(r) 
n(r) 



(r^ + a)- 1 



fi(r) (l-4a/r 3 / 2 + 3a 2 /r 2 ) V2 , 
n(r) (1 - 6/r + 8a/r 3/2 - 3a 2 /r 2 ) 1/2 



(2-1) 



the rotational, vertical epicyclic, and radial epicyclic frequencies, respectively. The angular 
momentum parameter a is less than unity in absolute value. 

The inner edge of the disk is at approximately the radius of the last stable free-particle 
circular orbit r = Tj(a), where the epicyclic frequency /«(rj) = 0. This radius r^(a) is a 
decreasing function of a, from r$(— 1) = 9 through rj(0) = 6 to r»(l) = 1. So all the relations 
we use are for r > rj, where «(r) > 0. Note, in particular, that 



f2(r) > Qj_(r) > n(r) , a > ; fi±(r) > f2(r) > «(r) , a < . 



(2-2) 



For a = (a non-rotating black hole), Q(r) = Q±(r) > k(t). We will also employ the 
outer disk radius, r Q . For stellar mass black holes in low-mass X-ray binaries, orbital sepa- 
rations require that r < 1O 5 M /M. The properties of the outer disk have not been greatly 
constrained from observations, but the matter impacting from the companion star is likely 
to significantly disturb the outer region of the disk. We neglect such complications to our 
equilibrium model in this exploratory treatment of the outer modes. 

To simplify the analysis, we here consider barotropic disks [p = p(p), vanishing buoyancy 

freque ncy; a generalizati on to a small non-ze ro buoyancy may be found in lSilbergleit. Wagoner fc Ortega-Rc 
(120011 ). section 4.3, and iPerez et al.l (Il997l )]. In this case hydrostatic equilibrium provides 
the vertical density and pressure profiles 



P = Po(r)(l - y 2 ) 9 , p = po(r)(l -y 2 ) 



2\S+1 



<7 = i/(r-i)>o 



(2-3) 



-4- 



where T > 1 is the adiabatic index (for brevity and convenience, we will use parameters g 
and T alternatively). One has T = 4/3 within any radiation pressure dominated region of 
the disk, and T = 5/3 within any gas pressure dominated region. The disk surfaces are at 
y — ±1, with y related to the vertical coordinate z by 



y 



h(r) 



r- i 

2T 



and h(r) specified by equation (2.19) of iPerez et al.l (119971 ). More information on the unper- 
turbed disk is given below. 



2.2. Equations for the disk perturbations 



To investigate the eigenmodes o f the d isk oscillations, we apply the general relativis- 
tic formalism that llpser fc Lindbloml (119921 ) developed for perturbations of purely rotating 
perfect fluids. Neglecting the self-gravity of the disk is usually a very good approximation, 
as will be discussed in section 8. Viscosity corrections to the perturbations are of order a* 
(the usual viscosity parameter, not to be confused with the function a introduced below) 
and are assumed small. (However, they are fully included in the equilibrium model via a*.) 
One can then express the Eulerian perturbations of all physical quantities through a single 
function 5V oc 8p/ p which satisfies a second-order partial differential equation. Due to the 
stationary and axisymmetric background, the angular and time dependences are factored 
out as 5V = V(r, z) exp[i(m0 + at)], where a is the eigenfrequency. Then the assumption of 
strong variation of modes in the radial direction (characteristic radial wavelength X r <C r) 
ensures the approximate WKB separability of variables in the partial differential equation 
for the functional amplitude V(r, z) = V r (r)V y (r,y). The function V y varies slowly with r. 

The resultin g ordinary different i al equations for the verti cal ( V v ) and radial (V r ) eigen- 

functi ons are [see lNowak fc Wagonerl (119921 ) . IPerez et al.l (119971 ) and lSilbergleit. Wagoner fc Ortega-Rodrigi 
d200lh for details]: 



9x CPVy (IVy n 



d 2 V r 
dr 2 



(iO 2 - K 2 ) 



d_ 

dr 



dV r 

dr 



> (1-73) 



+ a 2 (u; 2 



1 - 



Vy = 0, 

V r = . 



(2-4) 
(2-5) 



The coefficient a(r) is defined as 



a(r) = PreiVg^/c s (r,0) , (2-6) 
where f3 Te i = dt/dr and Kerr metric coefficient in Boyer-Lindquist coordinates. 
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As in ISilbergleit. Wagoner fc Ortega- Rodriguezl (120011 ) , we assume that 



a(r) = 7(7")- 



r - r, ; 



< n < 1/2 , v > 



(2-7) 



where 7(7*) is some function bounded from above and away from zero, varying slowly with 
radius, and tending at infinity to a limit 7(00) = 700 > 0. Thus for very large radii a(r) is 
effectively a nonnegative power, 

a(r) w 7oor" , (2-8) 

while near the inner edge of the disk 



a(r) ~ ati(r — r^)" 



«i = 7( r iJ?" 



;2-9) 



The singularity at r = is due to the fact that the pressure at the inner edge vanishes 
in a typical disk model [p oc (r — rj) fc ] if the torque does, in which case \i — (T — l)k/2T. 
Typically, k = 2 and T = 5/3, giving fj, — 2/5. If even a small torque is applied to the disk at 
the inner edge, which is likely due to magnetic stress (IHawley &: Krolikll200ll ). the pressure 
and the speed of sound become nonzero at r = r^. This corresponds to a nonsingular a(r) 
with fi = 0. 

It is also true that v < 1 in most disk models; if it is less than one-half, some limitation 
on the disk size (given 7^, and vice versa) applies. Indeed, as stated in section 2.1, 

h(r) ~ c s /Q ~ 1/aVL < r . 

This relation, according to formula f)2-8p for a large outer radius r Q ^> r^, converts into 



<7c 



(2-10) 



and we have the following three cases: 

1. For v > 1/2, the outer radius of the disk should satisfy r D ^> / y o^ u ■ 

2. For v — 1/2, it is not the size of the disk, but the speed of sound that is limited by 
the requirement 7^ 3> 1. 

3. For < v < 1/2 the outer radius is limited by r Q <ti 7^ ( ' ' 5 . 

One should specifically bear in mind this last limitation in the following analysis, since 
0<z/<l/2in most accretion disk models. 



To gether with the appropriate homogeneous bounda ry conditions [discussed by lPerez et al. 



( 119971 ). ISilbergleit. Wagoner fc Ortega- Rodriguezl (1200 ll ). and below in sections 3, 4], equa- 
tions ()2-4p and (12-5 j) generate the vertical and radial eigenvalue problems, respectively. The 
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radial boundary conditions depend on the type of mode and its capture zone, as discussed in 
the following sections. The coefficient a, the vertical eigenfunction V y and eigenvalue (sep- 
aration function) \1/ vary slowly with radius (comparable to that of the equilibrium model), 
as does the (dimensionless) ratio of the corotation frequency ou(r) and Q±(r): 



(Note that \I/ and a;* depend also on the eigenfrequency a.) 



(2-11) 



Our goal is to determine the vertical eigenvalues \l/(r, a) and the corresponding spectrum 
of eigenfrequencies a from the two eigenvalue problems for equations (12-4R and (12-51) . Along 
with the angular mode number m, we employ j and n for the vertical and radial mode 



numb ers (n umber of nodes in the corresponding eigenfunction ), respectively, as in lPerez et al. 
(119971 ) and ISilbergleit. Wagoner fc Ortega- Rodriguez! (j200ll ) . 



2.3. Classification of Modes 

One can see from the radial equation (12-51) that a mode oscillates in a domain of r > 
where the quantity (uj 2 — /t 2 )(l — $ /uj 2 ) is positive. The first factor here (for a in the 
allowed range) changes sign twice at the points r = r±(m, a, a), r% < r_ < r + , the Lindblad 
resonances determined by equations (12-11) and (12-lip. This factor is negative between r_ an d 
r + and positive otherwise [see fig. 1 in ISilbergleit. Wagoner fe Ortega- Rodriguez! (120011 )]. 



Thus the classification of modes is based on the behavior of the second factor, containing 
the eigenvalue. 

There are essentially three possibilities, corresponding to the following types of modes: 

• g-modes: >(typically 3>)<^; capture zone r_ < r < r + . 

• p-modes: $ < (typically <C) oj 2 \ capture zone rj < r < r_ or r + < r < r Q . 

• c-modes: \I/ = oj 2 , for at least one value of the radius. 

The capture, or trapping, zone is where the mode oscillations occur. 

The g (inertial-gravity)-modes are centered on the radius r m where r_ = r + , correspond- 
ing to the maximum eigenfrequency |er|. The corot ating eigenfrequen cies \u\ corresponding 



to the lowest mode numbers are close to ft(r m ) [see iPerez et al.l (119971 ) for details]. 



The c(corrugation)-modes which have so far been investigated are of a special sort: 
^ — u) 2 fa throughout their whole capture domain. They are non-radial (m ^ 0, typi- 
cally m = ±1) vertically incompressible waves near the inner edge of the corotating disk 
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that precess around the angular momentum of the bl ack hole. Their funda mental fre- 



quency is shown to coincide with the Lense-Thirring (ILense fc Thirrind Il918l ) frequency 
(produced by the dragging of inertial frames generated by the angular momentum of the 
black hole) at their outer capture zone boundary r r ( m, a) in the appropriate slow-rotation 
limit (jSilbergleit. Wagoner fc Ortega-Rodriguezll200ll ). 



In this paper we invest igate p (inertia l -press ure)-modes, which have been studied since 
the pioneering a nalysis of iKato fc Fukud (jl980l ). However, the only previous relativistic 
treatment was by iPerezI (119931 ). As mentioned in the introduction, in this paper we concen- 
trate on the properties of the fundamental (m = j = 0) low frequency p-modes. 



3. Perturbative Solution of the Vertical Eigenvalue Problem 

We will use the self-adjoint form of the vertical equation (l2-4p . 



d 

dy 



+ 2gu% 



1 



1 - 



out 



v y = o 



(3-1) 



For the lowest angular mode number m = we have the corotation frequency coinciding 
with the eigenfrequency, uj 2 = a 2 [from equation fl2-lip ]. Hence we can rewrite the above 
equation as 



CV y = - (* 2 /nl) 2g [1 - (1 - *M 2 ) (1 - y 2 )] (1 - y 2 ) 9 - 1 V y 



Here the Hermitian operator 



d 
dy 



2\9 



(i-9 2 ) 



dy 



(3-2) 



(3-3) 



acts on those smooth enough functions of y on the interval (—1, 1) which are bounded 
together with the ir derivative at y = ±1. The last requirement is dictated by physics 
(jPerez et al.l 119971 ). It also turns out to be a natural boundary condition for the singular 
ODE fl£Tj) . 

We are interested in this first paper in the low-frequency (LF) p-modes. We specify the 
LF range by requiring that the following frequency ratio be small enough throughout the 
appropriate capture zone: 

a 2 /nl(r,a)<l. (3-4) 

The frequency Q±(r,a) decreases monotonically with radius for the rotation parameter a 
within the interval — 1 < a < a± w 0.953; for a > a± it has a single maximum near the 
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inner edge, but is still decreasing far out as r~ 3 / 2 . Thus for the outer LF p-modes captured 
within r + < r < r D , the requirement (13-41) needs to be checked only at the outer edge of the 
disk (r = r D ). The LF inner fundamental p-modes are concentrated very close to the inner 
edge, so the requirement ()3-4p needs to be checked only at r = for them. 

Note that from equation (12- ip . 

n±(r i ,a)| a=1 =fi ± (l,l) = 0. 

Hence for very rapidly rotating black holes (a > a±) the inner modes we are studying are 
not well defined. Some more specific constraints on the parameters for both the outer and 
inner LF p-modes are given later in this and following sections. 

Having set up the basic small parameter (13-41) . we solve the vertical eigenvalue problem 
by means of perturbations, 

iT, / \ 2 / \ 4 



K »^ (0, + (si) V(1) + fe) V<2) + "" (3 " 6) 

with the coefficients Xk and functions V^(y), k = 0, 1, 2, . . . , to be determined. (Recall 
that both quantities are, in general, slowly varying functions of the radius, as is fij_.) The 
substitution of these expansions in equation (13-21) leads to the following chain of equations 

CV^ = , (3-7) 

CV^ = -f k ; (3-8) 



fk(y) = 2g 



k-l 



(i - yy- 1 v^\ y ) - (i - y y x^ k - l - l \ y ) 



k > 1 



with the boundary conditions that all V^ k \y), k > 0, and their first derivatives be finite 
at y 2 = 1. Note that the right hand side fk(y) is a linear combination of V^\y) whose 
coefficients depend on Xi, with < i < k — 1. 

The only solution to the described 'ground state' (a = 0) problem is 

V (0) {y) = constant = 1 . (3-9) 

This shows that, at least to lowest order, the vertical eigenfunction has no nodes (j = 0) and 
is even. A straightforward symmetry and induction argument immediately demonstrates 
that all the corrections V^ k \y), and hence the vertical eigenfunction V y (y), are also even. 
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This allows us to work on the halved interval y G (0, 1) with the proper symmetry condition 

i 

at y — 0. It is convenient to introduce the notation (/, g) = f f{y)g(y)dy for the scalar 

o 

product. 

The procedure of successive solutions of equations (13-81) is standard. For a given k > 1, 
provided that all the previous equations have been solved and the eigenvalue corrections 
Xo, Xii ■ ■ ■ iXk-2 already found, we first determine Xk-i from the solvability criterion (the 
operator L is Hermitian), 

(V<®,f k ) = -(V {0 \CV {k) ) = -{CV {0 \V {k) } = . (3-10) 

We then find V^ k \y) by direct integration of equation (!3-8p . with a now completely de- 
termined right hand side. Of course, V^ k \y) is originally defined up to an additive term 
proportional to V^(y) = const; it is usually made unique by means of the orthogonality 
condition 

(y(°),y(*)) = 0, k = l,2,... . (3-11) 

We begin by solving the problem for k = 1. According to expression (13-91) . the r.h.s. of 
equation (13-81) for the first correction is 

h{y) = 2g [l-Xo(l-Z/ 2 )] (l ~ , 
so condition (I3-I0p in this case becomes 

o = (y(°), fx) = 2 9 f [i - xo (i - y 2 )} (i - y 2 ) 9 ' 1 dy , (3-12) 

which allows us to easily calculate 

Xo = 1 + l/{2g) = (r + l)/2 . (3-13) 
^From the definition (13-31) of C, equation (13-81) for V^'(y) after one integration reduces to 



dy 



29 (i-y 2 y 9 Ai-xo (i-^ 2 )] (i-vY^dv 

Jy 



which is finite at y = 1. Also, due to equation (13-121) . the derivative goes to zero at y — 0, 
as it should by symmetry. It turns out that the last equation simplifies remarkably to just 
dV^ /dy = y, and the first correction to the eigenfunction satisfying condition (13- lip is thus 

VW(y) = 0.5(y 2 -l/3) . (3-14) 
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To see how the expansion (I3-5P is generated, we now want to calculate Xi] th & t ls i to 
explore the solvability criterion (I3-I0p for k = 2. According to equations (13-81) and (I3-I4p . 
the r.h.s. of the equation for V^(y) becomes 



hiy) =9[(y 2 - 1/3) [l - xo (i - y 2 )} (1 - y 2 ) 9 -2(1- y*) 9 Xl 

Using this and formula (I3-I3p . from equation (I3-I0p with k = 2 it is not very difficult to find 
that 

(T 



Xl 2g(2g + 3) 2 (3T - 1) ' (3 " 15) 

So the equalities ()3-5p . (13-61) and (13-131) — (13-151) provide together the two term expansions 
for both the vertical eigenvalue and eigenfunction of the LF fundamental p-modes, namely: 

2 ,„ 1)2 



I"' 



r + i 



+ 



(T 



2(3r-i; 
> 2 1 / 1 

(7 \ 1 / 2 1 



+ . . . 



(3-16) 



(3-17) 



Apart from everything else, the first of these expressions gives a plausible criterion of how 
small our expansion parameter (13-41) should be. For our applications it is sufficient to have 
the second term in the expansion (13-161) smaller than the leading one. That is, 



a 



xo = (r + i) (3r-i) 
xi (r-i) 2 



— > 15, i < r < 2 . 

Xi 



(3-18) 



As mentioned above, this should be verified for r = r Q in the case of outer modes and for 
r = Ti in the case of inner ones. 



So, to lowest order 



1 - 



r + i 



(3-19) 



uj; 2 

is independent of the radius (if r is). The y alue of this constant for T = 4/3 agrees with 
the previous numerical results of iPerezI (1l993l ). which we have checked by our own numerical 
integration. Note also that the vertical eigenvalue ^ defined by equation (13-161) is negative, 



r - i 



a 



(r-ir 
2 (3r - 1 



<0; 



and varies slowly with radius, since w*(r) and fl±(r) do so. 
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4. Transformation of the Radial Equation and Boundary Conditions 



It turns out that the radial equat i on (12-5 j) is not convenient for analytical solution, be 
itaWKB approach as in iPerez et al.l (119971 ) and ISilbergleit. Wagoner fc Ortega- Rodriguez 
(120011 ). or some other asymptotic methods used in the sequel. To convert it, we first recall 
its self-adjoint form, namely 



d 
dr 



dV r 



(uj 2 — k 2 ) dr 



+ a 2 1- 



V r = 



(4-1) 



and then [as in lPerez et al.l (119971 ). section 3.1] invoke a definition of a new eigenfunction 

1 dV r 



W 

In terms of W(r), equation (14- ip becomes 

V r : 



k?) dr 



(4-2) 



1 



dW 



a 2 (l-^/cu 2 ) dr ' 

and by substituting equation (I4-3P into the definition (14-21) we obtain the equation 

1 dW~ 



(4-3) 



d 
dr 



ct 2 (l-^/cj2) dr 



+ k - k )W = 



(4-4) 



for the new radial function. This is the exact equation for W(r), equivale nt to equatio n ( |4^TI) 
for V r (r). In contrast, an approximate equation for W(r) was used in IPerez et al.l (119971 ) 
[see formula (3-3) there]. It was obtained by neglecting the derivative of a 2 (1 — ^! /uj 2 ). 
Fortunately, both equations lead to the same result for g-modes to lowest WKB order. 

In the present case of p-modes we introduce the radial independent variable 

dr , r* = const . (4-5) 

)} 

Since by definition of the p-mode 1 — ^/uj 2 > in its capture zone, r(r) is a monotonic 
function of the radius there, so the inverse map from r to r is well defined. In terms of r, 
equation (14-41) reduces immediately to the standard WKB form that we need: 



rr 




/ «V) 









d 2 W 

~dp 2 



+ S(t)W = 



S(t) 



a- 



(4-6) 



[T his particula r reduc tion of a general second order ODE to the WKB form was suggested 
by iPonomarevi (119651 ).] The coefficient S(t) is positive in the inner, outer, and full-disk 
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p-mode capture zones, ri < r < r^, r + < r < r Q} and < r < r , respectively. The 
boundaries r± are the turning points of equation (14-61) . Expressed via the r variable, the 



corresponding capture zones are r» < r < r_ 



< r < r , and 7* < r < r . We use 



the natural notation r A = r(r A ) for A = ±, i, o. The function S(t) is completely specified, 
because a(r) oc c7 1 (r, 0), oj(r), w*(r) and «;(r) are given by equations (12-71) . (12- lip , and (12- ip . 
and r = r(r) is the inverse of the function (14-51) . 

We now move on to the radial boundary conditions. The physical processes and condi- 
tions at the edge of the disk are uncertain, but we assume that there is some reflectivity of 
the mode there so that a boundary condition is appropriate. We parameterize (as usual) the 
lack of knowledge of this boundary condition by setting the most general conditions there: 



a dVr -a t/ 
cos &i — sin &i V r 

dr 







a dVr -a t/ 
cos o — sin u Q V r 

dr 







with — 7r/2 < 9i >0 < 7r/2. Because of equations (14-21) . (14-31) and (14-51) . this translates into the 
following boundary conditions in terms of r and W{r): 



< 2 (ri) cos 9i W + sin Q { ^f- 

dr 







[uj 2 {r ) - K 2 {r )] cos O W + sin8 



dW 



dr 



(4-7) 



(4-8) 



recalling that n(ri) = 0. 



This is the complete set of boundary conditions for the radial equation (14-61) for the 
full-disk modes. However, for the inner and outer p-modes only one of these conditions is 
employed, since the other boundary is formed by one of the turning points, either r_ or r + . 
In both cases we require the solution to decay outside the turning point boundary of the 
trapping zone; that is, for r = r_ + in the inner and for r = r + — in the outer case. This 
requirement concludes the formulation of the eigenvalue problem for W(t) for p-modes of 
all types. 



5. Outer Low-Frequency Fundamental P— Modes 

Since m = and we assume a small eigenfrequency, r + « cr~ 2 / 3 . This turning point 
should be within the disk, so 



a 



-2/3 



) > n > l 



(5-1) 



13 



Thus, from equations (12- ip and (12-81) . for all points in the capture zone r + < r < r Q we have 
(to lowest order) 

k(v) « £l_i_(r) « fi(r) « r~ 3/2 , «(r) » 7 OCJ r i/ . (5-2) 
The inequalities (15-11) provide the eigenfrequency range 

r- 3 / 2 < a « rr 3/2 . 



However, the validity of condition (13-181) for the vertical eigenvalue approximation (13-191) at 
the outer edge gives, in view of equation (l5-2p . a much more restrictive upper bound on the 
eigenfrequency: 

1 



3/2 
To 



< a < 



y/(T + 1) (3r - 1) 1 



(r-i; 



3/2 



(5-3) 



For example, if T = 5/3 this reduces to 

r~ 3/2 < a < 5r~ 3/2 . 
Therefore the only situation which makes our assumptions consistent is 
a = r " 3 / 2 (l + e) , r + = r Q [l - 2e/3 + (9(e 2 )] , 



(5-4) 



with some small enough positive e to be found. Note that the upper bound on the eigenfre- 
quency in expression (15-3 j) stems from the use of the simplest approximation (13-181) . which 
nevertheless is valid if the eigenfrequency of the mode is low enough. Hence the whole fre- 
quency range (I5-3P is an immediate implication of the requirement that the eigenfrequency 
of the outer fundamental p-mode is low. 

For the outer p-modes the radial eigenvalue problem is considered on the interval r + < 
r < r . Hence it is convenient to chose r* = r + in equation (14-51) . so that 



T 



a 2 (r) 



r+ 



\]>(r ) 

ul{r') 



dr 



T+ = , T 



a 2 {r) 



r+ 



\I>(r 



dr 



(5-5) 



and the capture zone is < r < r Q . In our case (m = 0, low frequency) we replace the 
vertical eigenvalue with its expression (13-191) and a(r) with its outer edge value given by 
formula (12-81) to obtain 



2 2, 



2 2u I \ 

1oc r o (r-r+) , t 



r + i 



2 „l+2u. 



7oo r 



(5-6) 



(to lowest order in e). 
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We now need to solve the radial equation (I4-6P in the capture zone < r < t q with the 
boundary condition (14-81) at r = r G , and the requirement that the solution decays inwards 
through the disk for r < 0. By its definition (I4-6P and formula (15-61) . the coefficient in the 
radial equation simplifies to 



S(t) « s + t 



dS 
dr 



12 



r=0 (T + iy^r, 
Hence the radial equation (I4-6P becomes the Airy equation 

d 2 W 



4(l+i/) • 
o 



(5-7) 



dr 2 



+ S+tW = , < T < T , 



(5-8) 



whose solution satisfying the decay condition is the Airy function of the first kind, 

W oc Aii-s^r) . (5-9) 

Substitution of this eigenfunction into the boundary condition (14-81) . using equations (15-41) 
and ( 15-Th . provides the eigenfrequency equation in the form 



2e 



cos 6 Ai(— z) — sY 3 sin9 Ai (—z) = , z = sVV 



V3. 



"4(r + i)~ 


1/3 


Too 


9 




1/2-1/ 

lr J 



2/3 



(5-10) 

(A prime denotes the derivative in the whole argument.) Unless o is very close to zero, the 
first term here can be dropped, giving 



while in the opposite case 



Ai{-z) = , 
Ai(-z) = , 



e a = o. 



Both equations allow for an infinite sequence of positive roots < zo < Z\ < . . . < z n < 
which provide the eigenfrequencies according to equations ( 15-41) and ( 15-101) . namely: 



-n 









= r -3/2 (1 . 




"4(r + i)~ 


-1/3 


Too 


-2/3 


9 


9 




1/2-!/ 

Lrv 


.4(r + i)_ 



1 V3 



(5-11) 



Our requirement that e n be small enough (at least e n -C 4) can be easily met by its expression 
( 15- lip , since the last factor there is small due to inequality (12-101) . Since for large n the roots 
z n grow as 

Zn ss (3n7r/2) 2/3 , n>l, (5-12) 
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the largest possible radial mode number N is found from 



N < (2/3tt 
and can thus be rather large 



~4(r + i)~ 


1/2 


Too 


9 




1/2-u 

lr J 



c s {r o ,0)rl /2 Kr ) 



(5-13) 



Let us briefly summarize the physical properties of these outer LF fundamental p-modes. 
More quantitative information about them is given in section 8. 

1. The spectrum of the outer LF fundamental p-modes is finite, although the total 
number of them can be large. 

—3/2 

2. All the eigenfrequencies are only slightly greater than r , so that they are prac- 
tically independent of all other parameters, including the angular momentum parameter a. 
They have a rather weak dependence on the midplane speed of sound at the outer edge of 
the disk. 



3. The radial extent of these modes, r Q 



(2/3)e n r D , could be large compared 



to r,-. 



6. Inner Low-Frequency Fundamental P— Modes 

6.1. Radial eigenvalue problem 

For m = and low frequency a the outer boundary (r = r_) of the inner p-mode 
capture zone is close to the inner boundary (the inner edge r = of the disk). Thus it is 
natural to consider the eigenfrequency low if the linear approximation for the radial epicyclic 
frequency, 



n 2 {r) ~ q 2 {r — r^) 



q 2 = 



dn 2 
dr 



1/2 

ry -a 



3/2 / 3/2 

r„- I r ,• + a 



< 1 



(6-1) 



is valid throughout the whole trapping zone. In particular, at its outer boundary 

a 2 = K 2 {rJ) w g 2 (r_ — Tj) , r_ w j-j + cr 2 /q 2 . (6-2) 

Since the second derivative of K 2 (r) is of the order of the first (= q 2 ), the condition for the 
linear approximation (16- ip to hold is 



^7s 2 < i • 



(6-3) 
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That is our definition of the low frequency range for this case. In contrast with the inequality 
(16-31) . the condition (13-181) needed for the validity of the approximation (13-1 9[) for the vertical 
eigenvalue is easily met. It does not impose any serious constraints on the eigenfrequency, 
such as the lower limit (I5-3P in the case of the outer modes, at least when — 1 < a < a±_ « 
0.953. The function q(a) is plotted in Fig. 1. It increases monotonically within this whole 
range of the angular momentum parameter, from the minimum value q(— 1) ~ 0.014 up to 
the maximum q(a±) ~ 0.111. 



We now fix the variable r from equation (14-51) as 

dr ; n = 



rr 




/ «V) 






uJ 2 (r')_ 



T- 



a 2 {r) 



1 - 



qr(r ) 



dr 



(6-4) 



Then from equations (13-19}) . (12-91) . and (16-21) we obtain 



(r + 1) a? 
2(1 -2/i) 



n) 



1-2/x 



(6-5) 



Calculating the coefficient S(t) in the radial equation 
! representatio 

_|_ & 2 T 2/i/(l-2/Lt) 



2(1-2/1) \qj 

6j) using all the above approxima- 



tions as well as the representation (12-91) for a(r) near the inner edge, we obtain 
d 2 W 



where 



dr 2 
b 2 



T 



1/(1-2/.) 



T 



1/(1-2 M ) 



W = 0, < r < r_ 



(6-6) 
(6-7) 



feify) = (1 _ 2/i )(l+2M)/(l-2M) g 2 J 2/( p + 1)a 2]2/(l-2 M ) 

Equation (16-61) should be solved with the boundary condition (14- 7p and the condition of decay 
outside r_. Unfortunately, it cannot be integrated in terms of the known special functions, 
so we treat it asymptotically with respect to the parameter < fi < 1/2. 



6.2. Solution for a weak inner edge singularity (/i — > +0) 

We observe that for disks with a non-singular or weakly singular inner edge, to lowest 
order in /x <C 1/2 equation (16-6 j) becomes the Airy equation 



cPW 
~d^~ 



+ b 2 ( T _ - r) W = 0, = 6 2 (0) = [2g/(T + l)a, 2 ] 2 . 



Its solution satisfying the decay condition is W oc Az[— 6^ 3 (r_ — r)\. In complete similarity 
to the previous section, the boundary condition produces 

< 7- B = ^[2/(r + l)] 1/ V /8 «T 1/8 , n = 0,1,2,... , (6-8) 
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where z n is the (n + l)-th positive root of either Ai(— z) or Ai'(—z), depending on whether 
the function or its derivative has a larger coefficient in the boundary condition (14- 7p (i. e., on 
how cr 2 compares to b 2 J 3 ). The LF inequality (16-31) applied to the lowest mode from equation 
( 16-81) requires that 

(qa t ) 2/3 > 1 . (6-9) 

This is true for typical values of the parameters involved, so generically LF fundamental inner 
modes in the disks with a weak inner edge singularity do exist. In addition, the condition 
( 16-91) implies that cr 2 ^> b 2 J 3 \ therefore z n is the root of the Airy function Ai{— z) (unless 
Qi = tt/2). 

For large enough values of n the asymptotic behavior (15- 12[) for z n holds, while on the 
other hand the LF condition (16-31) must still be valid. Hence the largest possible radial mode 
number N is determined by 

N < (2g«,,/37r)[(r + l)/^ 1 / 2 , (6-10) 

for N>5. 

A physically meaningful way to represent the eigenfrequency is 

a n = z 3 J 2 [2/(r + l)] 1/2 c s (r_,0)// n , (6-11) 

where r_ = (r_) n , and l n = (r_) n — r« is the length of the mode capture zone. One derives 
this from equation (16-8 j) by removing q with the help of equation (16-21) . q 2 = a 2 /l n , and also 
noticing that for jii « and low frequency, the speed of sound c s (r, 0) ~ throughout 
most of the capture zone. 



6.3. Solution for a strong inner edge singularity 

The limit case of a disk with a strong inner edge singularity, 1/2 — 0, is apparently 
a singular one for equation (16-61) . To investigate it, we once again change the independent 
variable to 

x = (t/t^) 1 ^ 1 - 2 ^ (6-12) 

to obtain 

A ( x ^ d ^\ + ( 2 (l - x )w = , 0<x<l, (6-13) 
ax \ ax J 

with 
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To lowest order in 1 — 2/i — > +0, we can replace x 2fl with x in the equation (16-131) . Then we 
also make a substitution 

W(x) = e~ c - x w(x) , 

reducing it to 

x ^ + (1 - 2Cx) ^ - C(l - C) ™ = , < x < 1 . 

aar ax 

Up to a rescaling of the independent variable, this is the confluent hypergeometric equation, 
whose only solution non-singular at x = is the Kummer function (confluent hypergeometric 
function of the first kind). 

Hence 

(\-C \ (2Cr)^ 1+ ^/ 2 

W( X ) = e-f* 1, 2 C x) ~ ^ { >_ m *> - » , z - oo , 

unless (1 — C)/2 = — Ji, n = 0, 1, 2, ... [r_g(z) is the Euler gamma function]. The condition 
for the mode to decay on the right of x — 1 can be thus fulfilled in the latter case only. 
The Kummer function then becomes the Laguerre polynomial L n , and the eigenvalue and 
its eigenf unction become 

Cn = 2n + 1, W n (x) = e- {2n+1)x L n (2(2n + l)x) , n = 0,l,2,.... (6-15) 

The eigenfrequency is found from equations (I6-15P and (16-141) to be 

a n = {(2n + l) ^/(r+l)] 1 / 2 ? 2 ^)^} . (6-16) 



As in the previous case, for the inner LF fundamental p-modes in the disk with a strong 
singularity to exist at all, condition (I6-3P should hold at least with a from equation (I6-I6p . 
This implies that 

(g«,) 2/(3 - 2 ^ > 1 , (6-17) 

which is just slightly different from the corresponding inequality (16-91) . and typically is also 
valid. The maximum admissible mode number iV is determined by 

N < 0.5 [(r + 1)/2] 1/2 qon , (6-18) 

and its value can be large enough, comparable to the number of inner modes found for a 
weakly singular disk given by equation (I6-I0p . which is about two times smaller. 

Note that the expression (16-161) for the eigenfrequency to lowest order in (0.5 — fx) ~ 
simplifies to 

{(2n + l)[2/(r + l)] 1 / 2 (gM)} 1/2 . (6-19) 



Or, 
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In terms of physical parameters, that is the speed of sound and the trapping zone length, 
the eigenfrequency can also be written as 

a n = (2n + 1) [2/(r + 1)] 1/2 c s (r_, 0)/l n , (6-20) 

where, as before, r_ = (r_) n and l n = (r_) n — r». In this form it is very similar to the 
eigenfrequency (16-1 ip for the previous case, the only difference being the numerical factor 
(2n + 1) instead of zf/ 2 . 

In addition, if we formally set \x = in the expression (16-161) for the eigenfrequency 
[which, in fact, is valid in the opposite limit case (/x — > 0.5 — 0)], we obtain 

a n = (2n + 1)V 3 [2/{T + i)] 1 / 6 ^ 3 ^ 1 ' 3 , n = 0, 1, 2... . 

Quite remarkably, this expression essentially coincides with the 'small /i formula' (16-81) . with 
the exception of the first factor giving the dependence on the mode number [(2n + 1) 1//3 in 

1 /2 

place of z n ]. For the lowest mode [n — 0), even these factors differ by only about 5%. 

Finally, in both limit cases of the weak and strong edge singularity, it is not difficult 
(although hardly needed) to calculate higher order corrections using an expansion of the 
type 

u e = exp(elnu) = 1 + elnu + 0.5 (elnu) 2 + ... , e — > , 

with u = t, e = 2^/(1 — 2/i) for // w 0; and w = x, e = 1 — 2/i for /i « 1/2 — 0. The 
weak logarithmic singularity at u = does not prevent the use of standard perturbation 
techniques, as carried out in section 3. 



7. LF Fundamental P— modes vs. Fundamental C— modes 



It is both natural and useful to discuss the features of the LF fundamental p-modes by 

compa ring them to the LF fundamental c-modes found in lSilbergleit. Wagoner fc Ortega-Rodriguez 
( I200ll ). For brevity, in this section we speak about just p- or c-modes, always meaning only 
low frequency fundamental modes. 



7.1. Eigenfrequency spectrum 

a) Angular and vertical mode numbers. They are m = j = for p-modes and m 2 = j = 1 
for c-modes. 

b) Radial mode number. In both cases, the number of the radial modes (N) is finite. 
Therefore, the whole eigenfrequency spectrum is finite. 
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7.2. Spectrum dependence on parameters 

a) Dependence on the angular momentum. The p-modes exist within the range — 1 < 
a < o± ~ 0.953. The c-modes exist only for corotating disks, in the range < a < a < a±. 
The cutoff value Oq ~ 10~ 5 — 10~ 3 depends on the adiabatic index T, the behavior of the 
speed of sound through the whole capture zone (in particular, on both /i and u), and slightly 
on the boundary condition parameter 9^. 

For the outer p-modes the dependence of eigenfrequencies on a is very weak (through 
the value of the speed of sound at the outer edge only). For the inner p-modes it is essential, 
but for c-modes it is very sharp, especially near the cutoff value a n for the mode number n 
(with a n > a n -i > . . . > a± > a ). The total number of the modes, N, depends on all the 
parameters, but changes most strongly with a in the case of c-modes. 

b) Dependence on the boundary condition parameter 0j. The dependence of p-modes on 
9i is very weak. The dependence of c-modes on the boundary condition parameter is also 
not strong. 

c) Dependence on the speed of sound (i. e., on /i, u, ati, joo)- The p-modes depend 
essentially on the behavior of the speed of sound near either the inner or outer edge of the 
disk; that is, on \i and CKj (inner modes) or v and (outer modes). The c-modes depend 
on all the parameters, including a considerable dependence on v. 

d) Lower bound for the eigenfrequencies. The eigenfrequencies are bounded from below 
in all cases: \o~ n \ > [2/ (r + l)} 1 ^ 6 q 2 ^ 3 a i ^ 3 for inner p-modes, \a n \ > r 3 ^ 2 for outer p-modes, 
and |<x n | > 2aor~ 3 for c-modes. Note that the lower bound for c-modes depends on the 
angular momentum parameter and the outer radius of the disk (in fact, the bound is the 
Lense-Thirring frequency at the outer edge). For inner p-modes, the bound depends on the 
angular momentum parameter, the behavior of the speed of sound at the inner edge, and 
weakly on the adiabatic index T. For outer p-modes the only dependence is on the outer 
radius. 



7.3. Mode capture zones 



The capture zone of both inner and outer p-modes shrinks as the eigenfrequency ap- 
proaches its lower limit. The radial extent of the capture zone of the c-modes increases as 
their eigenfrequencies decrease, and can cover the whole disk. The difference in the behav- 
ior is mainly due to the diff erence in the angular modes numbers, m = and m = ±1, 
respectively. [See Figure 1 in lSilbergleit. Wagoner fc Ortega-Rodriguezl (120011 ).] 
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7.4. Oscillation forms 

To lowest order, for p-modes 

5p ip) /p oc 5V (p) oc /^ ) (r)exp[z ( r^t] , 

while for c-modes 

5p {c) /p oc 5V {c) oc 2^ c) (r)exp[i(±0 + 4 c) t)] . 

Therefore p-modes are almost purely radial, since the oscillations of particles in the vertical 
direction are suppressed, 

C oc [d(6V {p) )/dz] « . 



\£ z is the vertical component of the Lagrangian displacement from equilibrium (IPerez et al. 



19971 ).] The vertical oscillations of the particles in a c-mode are independent of z. Particles 



on the midplane do not oscillate in the radial direction, 

C oc [d(5V ic) )/dr] ^ at z = 0. 



8. Numerical Results and Discussion 

8.1. The unperturbed model 

We now apply the above analysis to models of black hole accretion disks. We start by 
stating our assumptions about the unperturbed model, in addition to what has already been 
said in section 2. The results discussed below were obtained by assuming that gas pressure 
dominates in the disk, so that T = 5/3, and that the opacity is due to (optically thick) 
electron scattering. These conditions hold, regardless of the luminosity, at least near the 
inner edge and for r ^> 7-j. [For example, for L/LEdd = 0.1, M/M Q = 10, and (the usual 
viscosity parameter) a* = 0.1, the inner region where gas pressure dominates extends from 
Ti to 1.15rj for a = 0, and from r» to 1.079rj for a = 0.953.] 

It is notable that the interior structure of the (zero buoyancy) accretion disk enters our 
formulation only through the constant T and the function a(r), defined by equation f)2-6p . 
We modeled a(r) in the way described by equation (12-71) . The relativistic factor 

PreiVg^ = (1 + ar" 3 / 2 )[(l - 37- 1 + 2ar- 3 / 2 )(l - 2r^ + a 2 ^ 2 )}- 1 / 2 (8-1) 

that appears in the definition of a(r) is of order unity except for values of a close to 1 and 
simultaneously r close to « 1. It is shown in Figure 2. 
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The values fo r the speed of sound are derived from the fully relativistic results of 



Novikov fc Thornd (119731 ): 



C s /C 

c s /c 
c s /c 



) 1/5 (a,M/M )- 1 / 1 %- 1 / 5 (a)r- 9 / 20 , r > r t , 

-1/10-9/20-1/5/ 



8.61 x 10~ 3 {L/L Edd 

8MxlO- 3 (L/L Edd ) l / 5 (a*M/M e) —r- - -~(a, 
x^- 1 /5( a ,r l )I?- 1 / 10 (a,r i )/ 1 / 5 (a)(r - r,) 2 / 5 , r^r i , f i>0, 
8.61 x lO- 3 (L/L S(W ) 1 / 5 («,M/M )- 1 / lo rr 9 / 2 %-V5 (a) 

xi3- 1/5 (a,r i )p- 1/10 (a,r l )F 1/5 (a)(l - 3) 1/5 , r = r i; /i = . 



(8-2) 
(8-3) 
(8-4) 



(As we mentioned above, the cases /i > and \i = correspond to the absence and presence, 
respectively, of torque at the inner edge of the disk.) In these expressions, rj (shown in 
Figure 3) is the 'efficiency factor' which relates the mass accretion rate to the luminosity: 
L = r]Mc 2 . (It is the energy lost by a unit mass as it accretes from far away to the inner 
edge.) The constant 3 is the fraction of the inflowing specific angular momentum at rj that 
is absorbed by the black hole (so there is a torque at rj unless 3 = 1). The functions B and V 
are relativistic correction factors defined in lNovikov fc Thornd (119731 ). which approach unity 
at large r. The function / (shown in Fi gure 4) is defined such that Q = f(a)(r — rj) 2 t o 
lowest order in r — r^, where Q [defined in lNovikov fe Thornd (119731 ): iPage &: Thornd (11974] )] 
is the factor responsible for making a(r) singular at r«. In the presence of a torque at r^, 
Q(n) = F{a){l- (3), where 



F(a) = (1 - 2r; 



2 -2\l/2/ 



2ar- 3/2 + aVr 2 )(l - 3^ + 2ar i ° /z ) 



-3/2s_i/2 



-5) 



This function is shown in Figure 5. 



one obtains v = 9/20 (which we employ for 



By comparing expressions (12-71) and 
all cases). Equation (I8-3P yields /x = 2/5 for the case of no torque at the inner edge, while 
equation (18-41) is consistent with /x = for the case of a nonvanishing torque at the inner 
edge. 

Moreover, we are now ready to determine the actual parameters involved in the eigen- 
frequencies, that is, ai and 700. According to the definitions (12-61) and (12-91) . 

^ = lim (r» - r)^/3 reiv ^/c s (r, 0) . 

r— *ri+0 

By employing equations (18-31) and (18-41) . we then obtain 

ai(a,M,L/L Edd ,a*) ■■ 



116(L/L i?(M )- 1 / 5 (a,M/M ) 1 / lo (/5 re/v ^;) ri . ? 
xS 1 /5 (a;ri) pi/io (a?rt)/ -i/5 (a)) ^ = 2/5 



(8-6) 
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and 

Oi(a, M, L/L Edd , a,) = 116(L/L i?dd )- 1 / 5 (a,M/M ) 1 / lo (/3 re/ ^)nr, 9/2 % 1/5 (a) 

xi3 1/5 (a,r l )P 1/10 (a,r i )F- 1/5 (a)(l - /?)~ 1/5 , pt = . (8-7) 

These two functions are shown in Figures 6 and 7. We see that in both cases «j has a value 
which is approximately constant except for a > 0.7. These properties are reflected in the 
properties of inner fundamental LF p-modes described below. 

Similarly, from equations (12-61) and (12-81) we obtain 

7oo = lim r~ u p rely /g^/c s (r,0) , 

r— >oo 

which, in view of equation (18-21) . gives [y = 9/20): 

7oo(a, M, L/L Edd , a,) = 116(L/ 'L EM )-V 5 (a*M / 'M ) 1 / lo (/5 re/v ^;) ri r / 1 / 5 (a) . (8-8) 



This function (which determines the properties of outer fundamental LF p-modes) is shown 
in Figure 8. 

Typical values for disk density and central object mass make the neglect of self-gravity 
a common assumption. An order-of-magnitude analysis shows that self-gravitational effects 
can be ignored, either at the unperturbed or the perturbed level (Cowling approximation), 
whenever 

r 10 11 , . 

<C (8-91 

GM/c 2 (M/M©) - 96 ' 1 J 

in normal units. These estimates assume a* = 0.1 and L = L Edd , with this condition (18-91) 
being only slightly less restrictive for lower values of the luminosity. 

Thus, for stellar mass black holes the condition is satisfied by many orders of magnitude, 
while one needs to be more careful with supermassive black holes when studying the outer 
modes. A black hole mass of 10 8 M Q , for example, implies the constraint r < 2x 10 3 GM/c 2 . 



8.2. Numerical results 

We present our major observationally relevant results in Figures 9-15 and Table 1, 
obtained from the formulas derived in sections 5 and 6. For the outer p-modes, we take 
r a = 10 4 (GM/c 2 ) for the stellar case and r Q = 10 3 (GM/c 2 ) for the supermassive case (unless 
otherwise indicated), and o ^ in the boundary condition. For the figures, we have chosen 
the radial mode number n = 0. As with the g- and c-modes, one might expect that the 
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lowest radial mode would be the one most easily excited and with the largest net modulation. 
The viscous parameter is taken to be a* = 0.1 throughout. 

In Figure 9 we plot the fractional radial extent of the fundamental outer p-mode as a 
function of the angular momentum of the black hole. This quantity is 2eo/3 [as defined in 
equation (15-111) ]. so we see that the model is self-consistent in that eo <C 4. The frequency is 

3/9 

then essentially just r . Note that given the dependence of the speed of sound on L and 
M [expressions (!8-2p . (18-31) and (18-41) ]. all the p-mode plots will show the same 'universal' 
characteristic: higher values of the frequency and fractional radial extent for higher values 
of L and lower values of M. 

Figure 10 shows the dependence of the eigenfrequency of the fundamental \i = 2/5 (no 
torque) inner p-mode on a. The fractional radial extent of this type of mode is shown in 
Figure 11. All of the modes satisfy condition (I6-3P . 

In Figure 12 we plot the eigenfrequ ency of the fundamental fi = inner p-mode as a 



function of a. Employing the results of lHawley fc Kroliki (120011 ). a value of j3 = 0.95 has 



been chosen. The fractional radial extent of this type of mode appears in Figure 13. Note 
that these (/i = 0) modes for which M = 1OM and L = LiEdd approach the limit of validity 
given by equations (16-21) and (16-31) . since ol/q 2 = 2.1(gaj)~ 2 / 3 = r_ — r; L — ► 0.66 as a — > — 1. 
Comparing Figures 10 and 12, we notice that the frequencies are less sensitive to the value 
of than is the radial extent. 

Figure 14 shows the total number of radial modes allowable by our approach [which 
assumes equation (16-31) . 

Another observationally relevant relation is that between the frequency and the lumi- 
nosity. For the case of the inner p-modes, such a relation can be obtained analytically: 

d\oga/d\ogL = 1/15 , (8-10) 

valid for both the /i — > and the /i — > 1/2 cases. 

Table 1 shows the frequency of the first three excited radial modes relative to the 
fundamental one. For the inner modes, this ratio does not depend on L, M, a, or any other 
parameter except f3, which we have taken to be 1.00 (corresponding to fi = 2/5) and 0.95 
(for n = 0). For the outer modes there is some dependence on the other parameters (the 
ranges of L and M are the same as in the figures). 
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8.3. Discussion 



The frequencies of all the fundamental diskoseismic modes that have been studied are 
collected in Figure 15. They have been compared with observations of pairs of 'stable' 
high frequenc y quasi-periodic oscillations in two low-mass X- ray binaries with black hole 
candidates by lWagoner. Silbergleit fc Ortega- Rodriguezl (120011 ). Identifying the oscillations 
with the g- and c-modes produced possible values of the mass and angular momentum of the 
black hole. The top curve in Figure 15 corresponds to the orbital frequency of a free 'blob' at 
the inner edge of the accretion disk. Only the (more realistic) /i = inner p-modes have been 
included in the plot. The outer p-mode curve represents an upper limit to the frequency of 
this type of mode. It corresponds to L — L Edd , M = 10M Q , and r Q = 10 3 GM/c 2 . 



One of the reasons that the (inner) p-mode was excluded as a source of the observed 
oscillations is its relatively narrow radial extent, as indicated in Figures 11 and 13. It lies 
very close to the inner edge of the disk, so has very little luminosity to modulate. In addition, 
the radial velocity within the equilibrium disk (which we have neglected) becomes significant 
near the inner edge. As indicated in section 1, we are now investigating higher frequency 
p-m odes, which cover more of the disk. We h ope to relate this work t o the traveling p-waves 
that iHonma. Matsumoto fc Katol (jl992j) and iMilsom fc Taaml (119971 ) claim to be generated 
in their numerical simulations of viscous accretion disks. 

The c-mode has a somewhat larger radial extent. The g-mode has the largest radial 
extent, and is the only high-frequency oscillation that does not extend to the inner edge of 
the disk. Given the uncertainties of the physical conditions there, it is the most robust mode 
(being trapped purely gravitation ally). The viscous damping and grow th of all three classes 
of modes has been investigated by lOrtega- Rodriguez fc Wagoner! (120001 ) . However, attention 
has not been focused on the outer p-modes because of the uncertain physical conditions near 
the outer edge of the disk. 

The inner p-mode could be an important source of modulation of the accretion flow. 
This could be most observable in the emission from the surface of (weakly magnetized) 
neutron stars in binary systems. In general, we might learn more about the uncertain 
conditions at the inner edge of the disk via these modes. 

In the future we hope to extend this perturbation analysis to more realistic models of 
accretion disks. I t prom ises to be challenging, since recent numerical simulations [such as by 
Hawley fc Kroliki (120011 )] indicate that a stationary equilibrium may not exist on the length 
scales of our modes. If so, one can still hope that its characteristic frequencies are less than 
those of our modes. We do plan to include the effects of radial velocity as well as greater 
disk thickness. 
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Table 1. Value of a n /a 
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Fig. 1. — Dependence of q, denned by q 2 = (dn 2 / dr) rv on the angular momentum of the 
black hole. 
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Fig. 4. — Dependence of / on the angular momentum of the black hole. It is defined in such 
a way that, for the case of no torque at Q = f(a)(r — r^) 2 to lowest order in r — Tj. 
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Fig. 5. — Dependence of F on the angular momentum of the black hole. It enters the relation 
Q( ri ) = F(a)(l-0). 
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Fig. 6. — Dependence of ctj on the angular momentum of the black hole, for the \x = 2/5 
case. 
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Fig. 9. — Dependence of the fractional radial extent of fundamental outer p-modes on the 
angular momentum of the black hole, for (A) M = 1OM , L = L Edd ,r Q = 10 4 GM/c 2 , (B) 
M = 1OM ,L = 0ML Edd ,r o = 10 4 GM/c\ (C) M = 1O 9 M ,L = L Edd ,r = 10 3 GM/c 2 , 
and (D) M = 1O 9 M ,L = 0.01L Edd ,r o = 10 3 GM/c 2 . 
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Fig. 10. — Dependence of the eigenfrequency of \i = 2/5 (no torque) fundamental inner 
p-modes on the angular momentum of the black hole, for (A) M = 10M Q ,L = Lem-, (B) 
M = 1OM , L = OMLem, (C) M = 1O 9 M , L = L Edd , and (D) M = 1O 9 M , L = 0ML Edd . 
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Fig. 11. — Dependence of the fractional radial extent of fi = 2/5 (no torque) fundamental 
inner p-modes on the angular momentum of the black hole, for (A) M = 10M Q , L = Lem-, 
(B) M = 1OM ,L = 0.01 W (C) M = 1O 9 M ,L = L Edd , and (D) M = 1O 9 M ,L = 
0ML Edd . 
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Fig. 12. — Dependence of the eigenfrequency of /i — 0, (3 — 0.95 fundamental inner p- 
modes on the angular momentum of the black hole, for (A) M = 10M Q ,L = Lem-, (B) 
M = 1OM , L = 0ML Edd , (C) M = 1O 9 M , L = L Edd , and (D) M = 1O 9 M , L = 0.01L £dd . 
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Fig. 13. — Dependence of the fractional radial extent of /j, — 0, (3 — 0.95 fundamental inner 
p-modes on the angular omentum of the black hole, for (A) M = 10M Q ,L = Lem-, (B) 
M = 1OM , L = 0.01W (C) M = 1O 9 M , L = L Edd , and (D) M = 1O 9 M , L = 0.01 W 
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Fig. 14. — Dependence of No on the angular momentum of the black hole; the number of 
radial modes iV (described by our model) satisfies iV <C N . The three curves correspond 
to the outer (broken line), \i = 2/5 inner (dotted line), and /i — 0,(3 — 0.95 inner p-modes 
(continuous line). A value of r Q = 10 4 GM/c 2 was used for the outer modes. 
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Fig. 15. — Dependence of characteristic disk fundamental frequencies on the angular mo- 
mentum of the black hole. The top one is the orbital frequency of a 'blob' at the inner edge 
of the accretion disk. The g- and c-mode curves bound the region corresponding to the 
parameters ranges 0.01 < L/LEdd < 1 and 10 < M/M Q < 10 9 . The two uppermost inner 
p-mode curves bound the region corresponding to 0.01 < L/L E dd < 1 and M = 10M Q , while 
the bottom ones bound the region corresponding to 0.01 < L/ L Edd < 1 and M = 1O 9 M . 
In all cases, /i = with (3 = 0.95. The outer p-mode curve corresponds to L — L E dd, 
M = 1OM , and r Q = WGM/c 2 . 



